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Exact results for the classical and quantum system of two 
vertically coupled two-dimensional single electron quantum 
dots are obtained as a function of the interatomic distance (d) 
and with perpendicular magnetic field. The classical system 
exhibits a second order structural transition as a function of 
d which is smeared out and shifted to lower d values in the 
quantum case. The spin-singlet <-> spin-triplet oscillations 
are shifted to larger magnetic fields with increasing d and are 
quenched for a sufficiently large interatomic distance. 

PACS numbers: 73.20.Dx, 36.40.Ei, 64.90.+b 



Quantum dots are artificial atoms which have been a 
subject of intense theoretical and experimental research 
in recent years Q. A very recent development is the 
study of vertically coupled quantum dots || or artificial 
molecules. In Ref. H a double dot system containing 
three spin-polarized electrons was investigated. A se- 
quence of angular momentum magic numbers was found 
which depends on the strength of the interdot tunneling. 
An analytically solvable model of an arbitrary number 
of vertically coupled quantum dots containing each two 
electrons but with a 1/r 2 electron-electron interaction 
and with neglect of the interdot tunneling was presented 
in Ref. . The first experimental realization of vertically 
coupled quantum dot structures were recently reported 

Two vertically coupled classical dots were investigated 
by two of the present authors as a function of the in- 
teratomic distance (d) [Q. We found first and second 
order transitions which are similar to structural phase 
transitions. For first (second) order transitions the first 
(second) derivative of the energy with respect to d is dis- 
continuous, the radial position of the particles changes 
discontinuously (continuously), and the frequency of the 
eigenmodes exhibits a jump (one mode becomes soft, i.e., 
its frequency becomes zero) . Here we will investigate the 
influence of quantum effects on the above structural tran- 
sitions. As an example we limit ourselves to the exact 
solvable artificial molecule with in each dot one electron. 
In the classical limit this system exhibits a second order 
transition. We also investigate how the spin-singlet <-> 
spin-triplet oscillations are modified by changing the in- 
teratomic distance. The equivalent two electron quantum 
dot, i.e. d = 0, was solved in Refs. We will show 

that the inter-layer distance introduces an additional pa- 
rameter with which the electron-electron interaction can 



be altered, and which modifies the behavior of the system 
considerably. 

The present system consists of two layers, separated in 
the ^-direction by a distance d, with in each layer [xy- 
plane) one electron which is confined by a parabolic po- 
tential. We take the same confinement strength in both 
layers and the barrier between the two layers is taken 
sufficiently large such that there is no tunneling between 
them. Within the effective mass approximation, this sys- 
tem is described by the following Hamiltonian 
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with m the effective electron mass, e the dielectric con- 
stant of the medium the electrons are moving in, and 
loq the radial confinement frequency. A magnetic field is 
applied perpendicular to the dots, B = Be Zl which in 
the symmetric gauge is described by the vector poten- 
tial A(r) = [B x f]/2, and the above Hamiltonian can be 
rewritten as follows 



H t = ^-(p 2 1 +p 2 2 )+U(r 1 ,r 2 ) 
2m 
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where the symbol u> c = eB/mc stands for the cyclotron 
frequency and the potential U is given by expression (|^) 
with the scaled confinement frequency ujq = 7^0 where 

7 = y/T+^U^I- 

A crucial property of parabolic confinement potentials 
is that the center-of-mass and the relative motion are 
separable. The center-of-mass coordinate and the rela- 
tive coordinate are R = (fi + r 2 )/2 and r = r\ — r 2 , 
respectively, with the conjugate momenta P = pi + p 2 
and p — (pi — p 2 )/2. The center-of-mass and relative 
motions are described by the following Hamiltonians 



W(P, R) = -^-P 2 + muj 2 R 2 - \uj c \P x R] z , 
4m 2 

H(p, r) = ^-p 2 + U(r)- ^uj c [p x r\ z , 



where 
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and U is obtained from U by replacing luq by Qq. 

In the following we do not have to consider the center- 
of-mass motion because it is a harmonic oscillator with 
ground-state energy and excitation frequencies which do 
not depend on the interatomic distance d. We concen- 
trate our attention on the relative motion. 

First we consider the classical problem. The Newton 
equation of motion for this system is 



r = \Ur . x r 



-VC/(r), 
m 
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with uj c — uj c e z . The classical ground-state configura- 
tion is obtained by minimizing the potential U(f). A 
second order transistion occurs when the Coulomb in- 
teraction energy and the confinement energy are of the 
same order. To emphasize this we write the Hamiltonian 
in dimensionless form where we express the coordinates 
and energy in the following units r' = (e 2 /e) 1 / 3 a -1 / 3 and 
E' = (e 2 /e) 2 / 3 a 1 / 3 , respectively, with a = mu> 2 /2. The 
potential reduces to 



U(r) 



(7) 



which after minimalisation leads to the equilibrium ra- 
dius for the ground-state configuration ro(d) = yl — d 2 
for d < 1, and ro(rf) = for d > 1. The corresponding 
energy, and its first and second derivatives with respect 
to the distance d arc E = (3 — d 2 )/2, dE/dd — —d and 
d 2 E/dd 2 = -1 for d < 1, and E = 1/d, dE/dd = -l/d 2 
and d 2 E/dd 2 = 2/d 3 for d > 1. The first and second 
derivatives are shown in Fig. [l] by the dashed curves. It 
is clear that there is a second order transition at d = 1. 
This second order transition is of the same type as the 
single second order transitions we found previously for 
the artificial molecules with twice six, seven and eight 
electrons (j. 

Next we calculate the frequencies of the eigenmodes. 
We express the frequencies in the unit u' = u>q/v2 and 
the equation of motion becomes 



uj r x r 



- 2V£/(r), 



(8) 



where the dimensionless potential (|7j) is used. To find the 
spectrum we take the Fourier transform of this Newton 
equation which leads to the following excitation frequen- 
cies 



uj r = 0, u>t = \J 6(1 — d 2 ) + 1 



d< 1: 



UJ± = y/2{l - l/d?) + (uj c /2) 2 ± w c /2, d>l. 



(9) 



The frequency u> r = corresponds to the rotation of 
the molecule as a whole, and u>b is the frequency of the 
breathing mode. The frequencies uj + and w_ correspond 
with the out-of-phase motions of the electrons around the 
center of the dots. These classical frequencies are shown 
in Fig. H for uj c = by the dashed curves. Note that 



one of the eigenmodes softens at d = 1 which induces the 
second order transition. 

Now we turn our attention to the quantum mechani- 
cal problem and investigate how quantum effects influ- 
ence the classical transition. In the quantum case the 
ground-state energy can be obtained as a solution of the 
Schrodinger equation. Due to the cylindrical symmetry 
of the problem it reduces to the following equation for 
the radial wave function (\&(r) = cxp(— imip)R m (r)) 
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and A = Hloq / \ar' 2 y/2) , which is the ratio of the quan- 
tum energy with the energy at the transition point. We 
solved this differential equation numerically, using a non 
uniform space grid as was done in Ref. |p"lj| . 

To see how quantum effects influence the classical tran- 
sition at d = 1 we solve Eq. ([h]) for ui c = and we take A 
slightly different from 0. In Fig. |l|(a) the first and second 
derivative of the ground-state energy (i.e. n = 0, m = 0) 
with respect to d is given for A = 10~ 4 , together with 
the classical result. Fig. 0(b) shows the same results for 
A = 0.2. One notices that quantum fluctuations smear 
out the classical transition which happens already for 
very small values of A. In order to reveal the influence of 
the quantization on the eigenfrequencies obtained classi- 
cally we define them in the quantum case as the smallest 
differences between the excited state energies of Eq. ( |l(i| ) 
devided by A. This is shown in Fig. for A = 10~ 2 . In 
agreement with the disappearance of the classical tran- 
sition there is no softening of a frequencymode in the 
quantum case. 

The effectcs of these quantum fluctuations can be un- 
derstood from a simple qualitative consideration. As we 
are interested in the behavior of the system close to the 
transition point we use the following expansion of the 
potential ( |TT| ) 



U(r) 
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where we neglect the constant shift and introduced the 
parameter Ad = (d — do)j 8 / 3 , which characterizes the 
deviation from the exact transistion point do = j~ 2 / 3 w 
1 — lo 2 /12loq. Thus in the quantum case the magnetic 
field shifts the transition point to lower d values. 

After the above expansion we arrive at the following 
radial equation which describes the quantum behavior of 
the system for small A values and close to the transition 
point 
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Note that there are two parameters d and A which control 
the quantum problem while in the classical case there is 
only one parameter d which controls the closeness to the 
transition point. The role of the other parameter A can 
be understood by considering the following limiting cases, 
where we take oj c — for simplicity. In the case when 
the parameter Ad is not too small one may neglect the 
quartic term in Eq. ( [l3| ) and we obtain the Schrodinger 
equation for a harmonic oscillator. This leads to the fol- 
lowing excitation frequency estimation u ~ E / 'A ~ \J Ad. 
Consequently, not too close to the transition point the ex- 
citation frequency demonstrates the classical soft mode 
behavior. In the opposite limiting case close to the tran- 
sition point (small Ad values) the quadratic term can 
be neglected. Scaling the variable r — > A 1 ^ 3 r it can 
be estimated that the excitation frequencies are of or- 
der lj ~ E/X ~ A 1 / 3 . Therefore, in the quantum case 
the soft mode frequency does not tend to zero, as is also 
apparent in Fig. p[ Comparing the frequency estimations 
in both Ad regions we find that in the region |Ad| < A 2 / 3 
quantization influences the classical transition point be- 
havior. For A = 1CP 4 , as in Fig. |l|(a), this region is about 
Ad w 10 -2 , which agrees with the numerical results of 
Fig 0(a). 

From Fig. [l] we also notice that in the quantum case 
the transition is shifted to smaller d values. As a measure 
for this shift we take the interatomic distance value at 
which dE/dd attains its minimum value for the ground- 
state energy (n = 0, m = 0). This interatomic distance 



(thick curve in Fig. 
(dashed curve in Fig. 



3) coincides almost with the d value 
J) at which the radial wavefunction 
(inset of Fig. ||) becomes maximal in the center of the 
dot. The shift of the transition to smaller d values can 
also be induced by applying an external magnetic field 
(see Fig. ||). The plotted results for lu c = 5 are also for 
the state (n = 0, m — 0), which is not the ground-state 
anymore, except for large d values. Thus in the quantum 
case the electrons sit in the center of both dots for smaller 
d values than in the classical case, which causes the shift 
to lower d values. 

It is known that in quantum dots |j) and in the D~ 
problem Jl^| oscillations between spin-singlet and spin- 
triplet ground-states occur as a function of the magnetic 
field strength. Here we want to investigate how these 
oscillations change when the distance between the dots 
increases and how they are altered or modified by the 
presence of the classical transition. It is obvious that 
when the interatomic distance d increases, a higher mag- 
netic field is needed to induce these transitions and that 
for sufficiently large d no transitions occur. 

In Fig. 0(a) these spin-singlet «-> spin-triplet oscilla- 



tions are shown for d = 0.5 and A = 0.1. In Fig. [|(b) 
we plot the phase diagram for the system with A = 0.1. 
When d increases, the transitions are shifted to higher 
magnetic fields. The classical transition is seen as fol- 
lows: for d > 1 we found no spin-singlet <-> spin-triplet 
transitions anymore, while for d < 1 oscillations do occur 
and this for all values of A. As in the classical system, 
the interlayer correlations are sufficiently strong to im- 
pose the one atom properties for d < 1, and for d > 1 
the system behaves like two decoupled quantum dots. 
The spin-singlet <-> spin-triplet oscillations in quantum 
dots are related to the Wigner crystallization state of 
the quantum dot Jl3|,|l4| what in our case corresponds to 
the transition at d — 1 and which explains why for d > 1 
those oscillations are absent. 
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FIG. 1. The first and second derivative (full curves) of the 
ground-state energy for uj c = 0, i.e. n = 0, m = 0, for the 
quantum artificial molecule as a function of d for (a) A = 1CP 4 
and (b) A = 0.2. The classical result (A = 0) is shown by the 
dashed curves. 

FIG. 2. The quantum mechanical frequencies (full curves) 
for A = 1CP 2 and the classical frequencies (dashed curves) as 
a function of d and for ui c = 0. 

FIG. 3. The dashed curve is the d value for wich dE/dd 
is minimal, which is a measure for the shift of the transition. 
The full curve is the minimum d value at which the radial 
wavefunction (inset of figure) attains a maximum at r = 0. 
The results are shown for lu c = and lo c — 5. 

FIG. 4. (a) The lowest energy levels for the two-electron 
double quantum dot problem (d — 0.5, A = 0.1). The singlet 
phase m = (S = 0) at low magnetic fields is followed by 
m = 1 (S — 1), m = 2 (S — 0), . . . phases as the magnetic 
field strength is increased, (b) Phase diagram for singlet and 
triplet ground-states for A = 0.1. Only the first eight transi- 
tions are shown. When d > 1 no transitions occur. 
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Figure 2 
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Figure 3 
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Figure 4 
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